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Abstract. 

We consider a directed walk model of a homopolymer (in two dimensions) which 
is self-interacting and can undergo a collapse transition, subject to an applied tensile 
force. We review and interpret all the results already in the literature concerning the 
case where this force is in the preferred direction of the walk. We consider the force 
extension curves at different temperatures as well as the critical-force temperature 
curve. We demonstrate that this model can be analysed rigorously for all key quantities 
of interest even when there may not be explicit expressions for these quantities 
available. We show which of the techniques available can be extended to the full 
model, where the force has components in the preferred direction and the direction 
perpendicular to this. Whilst the solution of the generating function is available, its 
analysis is far more complicated and not all the rigorous techniques are available. 
However, many results can be extracted including the location of the critical point 
which gives the general critical-force temperature curve. Lastly, we generalise the 
model to a three-dimensional analogue and show that several key properties can be 
analysed if the force is restricted to the plane of preferred directions. 
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1. Introduction 

The development of atomic force microscopy and optical tweezers has allowed 
experimentalists to micro- manipulate individual polymer molecules {e.g. Bemis et al. 
(1999), Haupt et al. (2002), Gunari et al. (2007)), and this has led to a considerable 
body of theoretical work describing the response of a polymer to an applied force {e.g. 
Halperin and Zhuhna (1991), Cooke and WiUiams (2003), Rosa et al. (2003)). Several 
situations have been investigated, including pulling a polymer off a surface at which 
it is adsorbed, pulling a polymer from a preferred solvent to a less preferred solvent, 
pulling a copolymer which is localized at an interface between two immiscible liquids, 
and pulling a collapsed polymer (in a poor solvent) to an extended form. In this paper 
we shall be concerned with the latter problem. 

The stress-strain curve of a linear polymer in a poor solvent, being pulled in an 
AFM experiment, has been measured by several groups (Haupt et al. 2002, Gunari et al. 
2007). The force-extension curve shows a characteristic plateau. For forces below this 
critical value the polymer will be in a collapsed state, while for forces above the critical 
value the polymer will be stretched. The plateau region would seem to indicate a first- 
order phase transition. Indeed Grassberger and Hsu (2002) have studied self-avoiding 
walks with nearest-neighbour attraction and an applied force at low temperatures (poor 
solvents): they predicted a first-order phase transition in three dimensions. On the 
other hand, they see no sign of a first-order transition in two dimensions. 

A well-studied exactly solved model of poor solvent polymers is the self-interacting 
partially directed self-avoiding walk model (IPDSAW). We begin our discussion by 
noting that a partially directed self- avoiding walk on the square lattice without self- 
interaction is intrinsically anisotropic with a preferred direction so that the polymer's 
size scales proportionally to its length, and one perpendicular to this so that the 
polymer's size in this direction scales sub-linearly. As we shall consider the square 
lattice with the polymer oriented one way, the preferred direction will be the horizontal 
direction and the other direction will be the vertical direction. 

In the absence of an applied force the critical point of this model was found by 
Binder et al. (1990), and is expected to model the polymer collapse transition (or 9- 
point). The exact solution of the generating function was found by Brak et al. (1992) and 
its singularity structure was rigorously elucidated. Prellberg et al (1993) used recurrence 
relations to generate very long series to estimate the exponents and the scaling function 
for the phase transition. A second-order phase transition similar to the ^^-point was 
foimd. The tricritical nature of this transition was described by Owczarek et al. (1993) 
based upon small parameter expansions and calculations of a related version of the 
model (see below). 

Various calculations were made by Owczarek et al. (1993) that are worth noting. 
Firstly a semi-continuous version of the model was solved explicitly, building on works 
by Zwanzig and Lauritzen (1968 and 1970) on a related model, and showed the same 
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tricritical nature (identical exponents) as the lattice model (from the small parameter 
expansion and the numerical work of Prellberg et al. (1993)). A later calculation by 
Prellberg (1995) on the asymptotics of the generating function of the lattice model of 
staircase polygons enumerated by perimeter and area implicitly demonstrated that the 
scaling function and the exponents were also the same as in the fully discrete model: 
the g-Bessel functions involved are the same in the two lattice models. This was made 
exphcit recently by Owczarek and Prellberg (2007). The low temperature scahng of the 
partition function was found by Owczarek (1993) in the semi-continuous model, and 
once again this was the same as found by Prellberg et al (1993) numerically for the 
fully discrete case. Wc argue below that the uniform asymptotic expansion given by 
Owczarek and Prellberg (2007) implies that they are indeed similar. 

Secondly, Owczarek et al. (1993) generalised the fully discrete model to include a 
parameter that counts the horizontal span of the walks. The generating function was 
found by generalising the approach of Brak et al. (1992) for the no-force case. This added 
parameter is equivalent to considering a force applied in the preferred direction. The 
semi-continuous model also contained this parameter. In the semi-continuous model 
it was clear that this parameter did not change the nature of the transition. In the 
discrete model it was also implicit that the collapse transition was unchanged by this 
parameter (and confirmed by the work of Prellberg (1995)), and so unchanged by such 
an applied force. The connection to an applied force was however not made explicit in 
the paper. Rosa et al. (2003) studied the model numerically with the connection made 
explicit, and added the consideration of the end-to-end distance scaling function for the 
discrete model, corroborating the unchanged tricritical nature of the transition when an 
applied force in the preferred direction is added. They also plotted the critical-force- 
against-temperature curve for the discrete and semi-continuous models. Again the work 
of Owczarek and Prellberg (2007) makes this conclusion explicit. 

Thirdly, in the appendix of Owczarek et al. (1993) a further generalisation of the 
model was considered and the full generating function of this further generalisation 
evaluated exactly. In this generalisation the parameter for vertical steps was replaced 
by two parameters: one for steps in the positive vertical direction and one for steps in 
the negative vertical direction. No analysis of this generating function was attempted. 
We note that a recent paper by Kumar and Giri (2007) considered pulling in both 
directions though their focus was on finite size effects calculated numerically rather 
than thermodynamic transitions and exact results. 

In this paper we make explicit the connection to applying a force in the vertical 
(non-preferred) direction of the generalisation discussed above. However we go further 
and solve for the generating function along the surface in the parameter space which 
should include the transition point. Hence we find an explicit expression for the phase 
transition point and so are able to plot exactly the critical force against temperature 
curve. We also observe that one of the tricritical exponents is unchanged when this non- 
preferred force is applied. This indicates, though does not prove, that the transition may 
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remain second-order even in the case of this type of force. This is a httle unexpected 
as the force must change the end-to-end scahng of the high temperature phase when 
pulhng in the vertical direction from sub-hnear to hnear, which it does not when apphed 
in the horizontal direction. 

We also derive functional equations for the generating function and show how the 
solution of this general class of problem can be streamlined with this approach. 

Before we explain our work on the vertical pulling problem we summarise all the 
known results on the horizontal pulling problem making explicit the results in terms 
of the applied force. Moreover, we apply the rigorous techniques of Brak et al. (1992), 
which had only been applied to the case of no applied force, to the horizontal pulling 
problem. In addition we discuss the behaviour of the force-extension curves based upon 
the exact results. Essentially we bring together all the known results and extend them 
as necessary for the case of horizontal pulling of a partially directed polymer in two 
dimensions. 

At the end of this paper we consider a three-dimensional analogue of the model 
and show it has similar behaviour to its two-dimensional counterparts. 

2. Model and definitions 

Consider the square lattice and a self-avoiding walk that has one end fixed at the origin 
on that lattice. Now restrict the configurations considered to self-avoiding walks such 
that starting at the origin only steps in the (1,0), (0,1) and (0,-1) directions are 
permitted: such a walk is known as a partially directed self-avoiding walk (PDSAW). 
For convenience, we consider walks that have their first step in the horizontal direction. 
Let the total number of steps in the walk be n. We label the vertices of the walk 
i — 0,1,2, ... ,n. Let the number of horizontal steps be and the number of vertical 
steps be riy. To define our model we will need finer definitions, so let us define Uy^ to be 
the number of (0, 1) steps (positive vertical steps) and ny_ to be the number of (0, —1) 
steps (negative vertical steps). If the walk starts at the origin let the position of the 
other end-point be {sx, Sy) so that the span in the horizontal direction is Sx and the span 
of the walk in the vertical direction is Sy. We therefore have 

n = rix + Uy 

= nx + ny++ny_ (2.1) 

and 

— "^x 

Sy = riy^ - ny_ . (2.2) 

An example configuration along with the associated variables of our model is 
illustrated in figure 1. 

To define our model we add various energies and hence Boltzmann weights to the 
walk. First, any two occupied sites of the walk not adjacent in the walk though adjacent 
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Figure 1. An example of a partially directed walk (the bold black path) of length 
n = 21 with Ux = 8, Uy^ = 8, and ny_ = 5 and having six nearest-neighbour 'contacts' 
(shown as intertwined (red) curves) so m = 6. The horizontal span is s^^ = 8 while the 
vertical span is Sj^ = 3. One end is fixed at the origin while forces are applied to the 
other end (horizontal fx and vertical fy). 



on the lattice are denoted nearest-neighbour contacts or contacts: see figure 1. An energy 
— J is added for each such contact. We define a Boltzmann weight u = e^"^ associated 
with these contacts, where /3 = l/ksT and T is the absolute temperature. Without 
loss of generality we shall take the units of energy to be such that J = 1 and therefore 
a; = e^, except when we discuss pulling polymers without any self-interaction (u; = 1) 
where we will have J = 0. 

An external horizontal force fx pulling at the other end of the walk adds a 
Boltzmann weight h'^'^ with h = e^^"^ . An external vertical force fy pulling at the 
other end of the walk adds a Boltzmann weight v^y with v = e^^^ . 

The partition function Znifx, fy,l3) of the model for walks of length n > 1, where 
for later mathematical convenience the first step of the walk is a horizontal step, is 

Znif.Jy,/3)= Yl /,^^('^V^{^)u;-{^) , (2.3) 

ip is PDSAW of length n 

where m{(f) is the number of nearest-neighbour contacts in the PDSAW, (p. The 
generating function G{z; h, v, u) is 

oo 

G{z; h, v,uj) = Y, Znif., fy, f3)z^ , (2.4) 

71=1 

so z can be considered as a fugacity for the steps of the walk and the generating function 
as a "generalised partition function". We shall denote the radius of convergence of 
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G{z; h,v,u}) as a function of z as Zc{h,v,u}). The mean values of n, m, and Sy, are 
given by 

dlogG dlogG 
in) = z^^, im) = , 

dlogG dlogG 

^^^^ " ~dfr' ^^"^ " ''~d^' ^ ^ ^ 

We shall define the number of n-edge partially directed walks with m nearest- 
neighbour contacts, horizontal span Sx and vertical span Sy as dn{sx, Sy,m) so that 

Zn{fxJy,f3)= Yl dn{sx,Sy,m)h'^v'yu;"' . (2.6) 

When V — 1, that is fy — 0, we let 

hn{sx-, m)^Y dn{sx, Sy, m) , (2.7) 

Sy 

and so bn{sx,m) is the number of n-edge partially directed walks with m nearest- 
neighbour contacts and horizontal span s^- So 

Znif., 0,P) = Y^ b^{sx, m)}f-uj"' . (2.8) 

Sx,m. 

Let the number of partially directed walks of length n be 6„ so that 

6n = 2^ dn{sx-,Sy,m) = Ybn{sx-,m) . (2.9) 

It is advantageous when working with the generating function to define different 
variables. Let the generating function G{x,y-^,y^,u!) be defined as 

oo 

G{x, y+, y-,u) = J] x'^^^^V+^'+^^^y-"'"^^^^"'^^^ , (2-10) 

n=l if 

where the sum over <^ is over all PDSAW's of length n. Then making the substitutions 
X — hz 
y+ = vz 

y_ = z/v (2.11) 
demonstrates that 

G{hz, zv, z/v, u) — G{z; h, v, u). (2-12) 
Hence we have 

Zn{fxJy:P) = [z"']G{hz, ZV , Z / V ,Uj) 

It will be useful to define 

g+ = y^u) and g_ = y^ui (2-14) 
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and 

q^y/^, (2.15) 
and on making the substitutions (2.11) 

q = u}z . (2.16) 
We define the reduced hmiting free energy as 

<f., fy, P) = lim - log [Z„(/,, fy, /3)] (2.17) 

n—>-oo n 

where the existence of the hmit can be estabhshed by concatention arguments. The 
radius of convergence Zc(h,v,uj) can be related to the free energy K{fx,fy,P) as 

i^{fx,fy,f3) = -'^ogZc{h,v,uj). (2.18) 

It will turn out that there is a single phase transition where the free energy is 
singular as a function of f3. We shall denote the phase-transition inverse-temperature as 
/5* = f3*{fx, fy)- In general a superscript of t implies the critical value of a parameter, e.g. 
the critical horizontal force at a fixed temperature and zero vertical force is /* = /^(O, /3). 
However, we use a subscript t to denote the critical values of the fugacities. The 
value of the radius of convergence of G{z; h,v,u!) as a function of z at cUf — e^* is 
zl{h,v) ^ Zc{h,v,ut{h,v)). 

When fy — we shall set y — — and so the generating function we need to 
consider for the case where there is no vertical stretching force is G{x,y,y,u!), noting 
that G{hz, z, z, cu) — G{z; h, l,u}), where 

oo 

G{z- h,l,u;) = J2 Znif., 0, (3)z^ = Yl ^)^'^^'"^" • (2-19) 

n=l nySxjm 

Note that for fy — we have 

q = yuj, (2.20) 

which agrees with (2.16) when the substitutions (2.11) are made since y = z when 

y = y+ = y- (that is, v = i). 

In subsequent calculations it is convenient to define the generating function 
Gr{x, y^,y_, u) for walks that have r vertical steps immediately after the first horizontal 
step, where — oo < r < oo. It will also be convenient to introduce the generating 
function G'^{x,y^,y.^,u!) — Gr{x,y+,y-,u!) for walks that have r > steps after the 
first horizontal step in the positive y direction and G~{x,y+,y^,u!) — G-r{x,y+,y-,u!) 
that have r > steps after the first horizontal step in the negative y direction. Clearly, 

^0 = ^0 • 

In the functional equation section (section 7) we need 

oo 

F{p)= J2 Gr{x,y,y,uj)p\ (2.21) 
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oo 

F+{p) = J2 Gn^, y+, y-, ^)P' (2-22) 

r=0 

and 

oo 

F-{p) = Y,G;{x,y+,y.,u)p^ . (2.23) 

r=l 

Note the asymmetry of the summation index in the final two definitions. Note that 

oo 

G{x, y, y, uo) = F{1) = ^ Gr{x, y, u) (2.24) 

r=— oo 

and 

oo 

G{x,y+,y.,u) = F+{l) + F~{l)= Gr{x,y+,y.,u) . (2.25) 

r=— oo 

3. Pulling in the preferred direction (fy — 0) 

In this section we will begin discussing the case of horizontal pulling force where fy — 0. 
3.1. No self-interactions (uj — 1) 

Let us begin by considering the case of no interactions and no force so that J = fx = 
fy — 0. It is easy to see that hm„_,oo rT^ logftn — log(l + -\/2). Let us consider the case 
of no vertical force so that fy — and v — 1. As discussed in the previous section we 
are interested in the generating function G{z; h, 1, cu) — G{hz, z, z, uj). As defined above 
this generating function converges when z < Zc{h,l,uj). Recall that the free energy 
K{fx,0,f3) is given by k = — log^^ (assuming that the limit defining the free energy 
exists, which will be proved in a later section). 

When UJ = 1 (which corresponds to turning off the vertex-vertex interaction and 
hence to good solvent conditions), G{z; 1, h. 1) satisfies the equation 

G{z; h, 1, 1) = hz{l + G{z; h, 1, 1)) + (1 + G{z; h, 1, 1)) , (3.1) 

1 — z 

so that 

l + (3.2) 

From this we can readily calculate the ratio (s^) / (n) and take the thermodynamic limit 
by letting z ^ Zc{h, 1,1) where 

Vh^ + 6h + l-l-h 
Zc{h,l,l)^ ^ ■ (3.3) 

In figure 2 is a plot of limz^;^^{sx) / {n) against (3 fx- Note that the stress-strain curve is 
qualitatively the same as that found experimentally for polystyrene in a good solvent 
(toluene) (Gunari et al. 2007, figure 1). 

For UJ ^ 1 the situation is more difficult and we derive some results about the 
generating function in the following sections. 
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Figure 2. The dependence of lim {s^) / {n) on (3 fx when lo = \. 

Z — >Zc 

3.2. Convexity and continuity for uj ^ 1 

In this section we establish the existence of the thermodynamic hmit for the canonical 
problem, and prove convexity and continuity. This implies continuity of the phase 
boundary. 

It will be convenient to consider a subset of the walks counted by (^(z; /i, 1, cj). 
Let a„(s2:,m) be the number of partially directed walks with n edges, m contacts and 
x-span equal to which satisfy the additional constraint that their last step is in the 
east direction. (That is, both the first and last steps are east steps.) These walks can 
be concatenated by identifying the last vertex of one walk with the first vertex of the 
other walk which yields the inequality 

Defining the partition function 

An{h, uj) = Y^ a„(s^, m)h'''uj"' , (3.5) 

Sx,m 

this implies the super-multiplicative inequality 

Ani+n2 

{h,uj) > A„,^{h,uj)An^{h,uj). (3.6) 

Since An{h, uj) < S^u^h^, it follows that n^-^ log An{h, uj) is bounded above for h,uj < oo. 
Hence the super-multiplicative inequality above implies the existence of the limit 

lim n-' log Anih, iu) = ~K{f„ 13) (3.7) 
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for h,ij < oo. 

We now recall the definition of the partition function 

Znif,, 0, /3) = 5] bn{m, s,)u^h'^ (3.8) 

and define Bn{h, l,a;) = Zn{fx,0, P)- Since An+i{h,u!) — hBn{h,u!) it follows that 

lim log S„(/i, u) = «(/,, (3) = «(/,, 0, (3) (3.9) 

n^oo 

for h,u! < oo. 

To prove that K{fx,P) is a convex function of (3 and we note that Holder's 
inequality implies that 



^m,Sa: / \Tn,Sx 



2 



so that 



and hence 



> J]6n(m,s,)(v^)"^(V^)'^ (3.10) 



Bn{hi,UJi)Bn{h2,UJ2) > S„(V<^1<^2, a/M?^)^ (3.11) 



n Mog5„(/ii,a;i) +n Mog 5„(/i2, ^2) . _i , p / /Tir / ^ 

^ > n logS„(V%/i2, V^i^2)-(3.12) 

This shows that n"^ log Bn{h,u;) is a convex function of P and Pfx. (It is convex as a 
surface, not just separately convex in each variable.) Since the limit of a sequence of 
convex functions (when it exists) is a convex function, K{fx, 0, /3) is a convex function of 
/3 and Pf^. Hence K.{fx, 0, (3) is continuous and is differentiable almost everywhere. Since 
K{fx,0,P) = —\ogZc{h,l,ij), the singularity surface z — Zc{h,l,u!) is also continuous 
and differentiable almost everywhere. 



3.3. The generating function for ui ^ 1 

While the calculation of the generating function for 7^ 1 appears in Owczarek et al. 
(1993) (see section 4 of that paper) wc feel that it is worth summarizing and presenting 
in a slightly different way. In order to consider the case 0; 7^ 1 it was necessary to 
generalize a method originally due to Temperley (1956) and used for the case h = 1 hy 
Brak et al. (1992). 

For convenience let the partial generating function Gri^; h, 1,00) in the case where 
V — 1 (no vertical force) be denoted as 

gr^2Griz]h,l,Lj) r>l, (3.13) 
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and go = Go{z; h, l,uj). Now the generating function required is G{z; h, l,uj) which is 
the sum over all r as 



G{z; h,l,uj) = Gr{z; h,l,uj) . 



(3.14) 



Noting that G-r{z; h, l,uj) = Gr{z] h, l,uj), we have the partial generating functions gr 
satisfying the relations 



go = hz + hz{go + gi + . . .) = hz{l + G{z] h, 1, u)) 



(3.15) 



and 



gr = hz'^^ 2 + J](l + uj'')gk + {I + uj') 9k \ . r > 1. (3.16) 

\ fc=0 k=r+l J 

See figure 3. These results imply that g^ satisfies the recurrence 
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k=0 k = k=r+l k=r + l 



Figure 3. Walks counted by gr, with r > 0, eitiier contain a single horizontal 
bond (and so are an J or an ] ) or can be constructed by appending an J or ] 
configuration of bonds to a walk counted by gk- Summing over all these possibilities 
gives equation (3.16). 



gr+i — (z + ujz)gr + uj'^hz^^'^{uj — \)gr + OJZ^gr-l = . 
We have q = ujz so that (3.17) becomes 

gr+i - {z + q)gr + q''hz{q - z)gr + qzgr-i = . 



(3.17) 
(3.18) 
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This can be solved with the Ansatz 

oo 

gr = y'J2Pm{q)q^'^ (3.19) 

m=0 

which is a solution provided that 

-X{z + q)+qz^O (3.20) 

and 

^-(^^ = iXqr^-q)iXq--zf -'^'^- ^^'^'^ 
Since po{q) ~ 1 this gives 

^mj^m^m^^ g^m^m(m+l)/2 

Pm[Q) = / \ I. \ T-rm / > u T- [6.22) 

nfe=i(v-9) nfc=i(^r-^) 

The quadratic equation (3.20) has the two solutions Xi = z and X2 = q = urz and 
the general solution for g^^, > is 

Qr = A^gi'^ + A^gf^ (3.23) 

where 

Arguments similar to those in Brak et al. (1992) show that A2 — and Ai can be 
determined by noting that 

go = ^A,g^'^ ^hz + hzG{z; h, 1, a;) (3.25) 

and 

g^^ Aig[^^ ^a + bG{z;h,l,u) (3.26) 

where a — hz'^{2 + hz — ujhz) and h — hz^{l + hz + cu — cuhz). These simultaneous 
equations give 

G{z] h, 1, uj) = —J-: — ^ , (3.27) 

bgj,'> - 2hzg'i> 

which can be written as 

, , 2hz^(uj - l)g^^^ 

(^ - 



n(0 - V < \r\^ ^ - q) q ^ " mr o^N 



z{l + hz + uj- ujhz)gj^^ ~ 2g[^^ 
The solution given by Owczarck et al. (1993) was written as 



(3.28) 



1 + G{x,y,y,u}) ^ —— \ , n vn (^-2^) 

2{n{y,yuj,xy{uj - 1)) - 1) + [1 - u;){l - x) 
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where 

njf ,N H{y,q,qt) 

^^^'^'^) = ^(^ ^^-^'^ 

and 



^(^-^-^) = E rl^Li^ • (3-31) 



The two forms can be seen to be the same when one notes that 
1 + G{z; h,l,u!) — 1 + G{hz, z, z, cu) 



2(n(z, zu, hz'^iu - 1)) - 1) + (1 - a;)(l - hz) 
uj-1 



(3.32) 



{1 + hz + uj — uhz) — 2H{z, zuj, hz'^iuj — 1)) 
and that 

niz.zu.hz^iu-l))^^^ . (3.33) 

The impUcations for the singularity diagram of this expression for G{z] /i, 1, a;), and 
the phase transitions, will be considered in the next section. 

3.4.- Solution on the special surface z = l/uj 

Before we proceed to discuss the singularity diagram let us consider the solution on the 
special surface defined by g = 1 since the solution described by (3.28) is singular when 
z — XjijJ though the generating function is finite for large enough lo. 
The recurrence for the partial generating functions becomes 

t/r+i + \hz{\ - 2;) - (1 + z)\gr + zgr-\ = , (3.34) 

and now the simpler Ansatz of gr — C^^ can be used to find 



l + GiM.M-sn ,,^,,T^:^^,^^ - (3.35) 
When /i = 1 we recover the previously calculated result of Brak et al. (1992). 

4. The phase diagram for a horizontal force 

In order to analyse the phase transition structure we need to analyse the form of the 
singularity diagram as a function of the uo and h parameters. We derive a functional 
equation for a slight generalization of the quantity g'i^^ We introduce a parameter t and 
define the function 

°^ h"''UJ^"'' ( — ll „m(m+5)/2^m 

g{t; q,h,u;)^l + J2 rr- n kin (4-1) 

nfe=i(i-r)(i-r/w) 
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Note that g^^ is equal to zg{tq; q, h, u) evaluated at t = 1 while gl^^ — g{l; q, h, cu). 
The function g{t; q, h,uj) satisfies the functional equation 

g{t; q, h, u) + g(qH; q, h, u;)/u; 

= (1 + l/cu + {l/cu - l)hqH/uj)g{qt] q, h, cu). (4.2) 

Defining 

H{t; q, h, a;) = g{t; q, h, (jj)/g{tq; q, h, u) (4.3) 

gives 

H{t; q, h, a;) = (1 + l/u + {l/u - l)hqH/uj) - ^^"^ , (4.4) 

H[qt]q,h,u;) 

which leads to a continued fraction representation for H{t] q, h,uj). Note that 
and therefore 

A/ 7 -, \ (u — l)H(q,h,u!) ,^ ^, 

1 + G{z]h,l,u!) = ^ — ^ \ ' . (4.6) 

(1 + + a; — Oi)hz)H(q, h^u) — 2 

The function H{q, h, cu) is singular on the hyperbola ^ = = 1 for all values 
of h and the only other singularities in G come from the poles corresponding to 
the denominator of G being zero, i.e. determined by the solution of the equation 
{1 + hz -\- uj — ujhz)H{q, h,u) — 2 = 0. This line of poles intersects the hyperbola 
at a point determined by solving the equation {1 + hz + uu — ujhz)H{l, h,uj) — 2 = 0, 
where ^(1, h,uj) is determined by the quadratic equation 









1 + - + 




-1 






UJ 



H{l,h,uj) + - = , (4.7) 

to 



which comes from (4.4) on setting q = t = 1. The variables lv and h are related through 
the equation 

ah'^{uj-l)-2u{u + l)h + uj^-u;^ = (4.8) 
which implies that the critical value of h is given by 

ft, = (4.9) 

UJ — 1 

Let us consider qc{h,l,uj) = Zc{h,l,uj)uJ. We have that < gc < 1 for 
1 < UJ < ujt{h,l). Wc have argued that there is a pole in the generating function 
which means that there is a zero of -D = 1/(1 + G). That is, the function 

^ [(1 + hq/uj + uj-hq)- 2H{q/uj, q, hq\uj - l)/uj^)] 

obeys 

D{q„uj,h) = 0. (4.11) 
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We also have analyticity of ?i for < g < 1, shown by the ratio test. When cu — 1 there 
is a simple pole in the generating function. We have 

Z,(,,l..).l^(l±^^ (4.12) 



and so that 

dD{q,l,h) 



< (4.13) 



dq 

at q = qc{h, 1, 1). The derivative is negative for all < g < 1. 

We now prove that there is a simple pole for all 1 < a; < u!t{h, 1). For this, we 
compute 

dD{q,uj,h) h 2 2/ ^^/ 2^ u t a\ 

^ = 7^-^ - 1)/^^ ) ■ 4.14) 

dq uj uj — Idq 

Observe that g, t) is related to the generating function Gsc{x, y, q) of staircase 
polygons enumerated by width, height, and area. Equations (4.6) and (4.9) of Prellberg 
and Brak (1995) imply that 

Gsc{x,y,q) = y[n{qy,q,qx) - 1] . (4.15) 

It follows that 7i{q/w,q,hq'^{u! — is a power-series in q with non-negative 

coefficients. Therefore, its derivative with respect to q is positive. Hence 

< (4.16) 

oq 

and therefore there is a simple pole in the generating function at all values of 1 < a; < 
ujt{h, 1). 

Now, using a version of the implicit function theorem that gives analyticity (see 
chapter 6 of Krantz and Parks 2002) it follows that qdh, l,a;) is an analytic function of 
both UJ and h. 

To determine the detailed shape of the part of the phase boundary determined by 
the poles of G we solve the equation {1 + hz + u — ujhz)H — 2 = numerically by 
evaluating H from its continued fraction expansion. That is, wc write H as 

H{t]UJ,h,q) = {l + l/uj+{l/uj-l)hqH/uj){l + ^oC) (4.17) 



where 



1 + 



and 



-1 
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To evaluate the continued fraction C efficiently we note that if Cm is the truncation 
of C at order m then Cm can be written as a rational function whose numerator and 
denominator depend on m and are determined by recurrences. This gives a convenient 
way to evaluate C to the required accuracy. In figure 4 we show the phase boundaries 
in the [uj, z)-p\a.ne for several values of h. 

0.5 n 



0.4- 



0.3 - 
0.2- 



0.1 - 



H • 1 ' 1 • 1 • 1 • 1 • 1 • 1 

1 2 3 4 5 6 7 

UJ 

Figure 4. The radius of convergence of the generating function as a function of lu 
when the force is applied in the x-direction. The rectangular hyperbola is independent 
of the value of h. The three curves with points marked correspond to /i = 1 (top 
curve), h = 1.5 and h = 2.0. 




5. Extension under a fixed horizontal force 

We shall now consider the average extension of the polymer as a function of the applied 
force at different temperatures. 

From section 3 we know that K{fx,0,f3) is a continuous function of (3 and fx- 
Moreover for any fixed fx the free energy ^{fx, 0, (3) is analytic for < /3 < oo except at 
at most one point which we label (3^ {fx) when it exists. This transition point is defined 
(see Owczarek and Prellberg 2007) via 

recalling that h = e^^'^. Similarly the free energy K{fx,0,f3) is an analytic function of 
fx except at the solutions of l3^{fx,0) = (3. Consideration of the equation (5.1) implies 
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there is a single solution. Let the solution for this 'critical force' be labelled /* and so 
by (5.1) it is given by the function 

which is the result (11) of Rosa et al. (2003). In figure 5 we plot the critical force /* 
against p'^. We note that for /5 < /3*(0,0) ^ 1.218 (equivalently T* ^ 0.8205) there is 
no positive critical force, and that for (3 > /5*(0, 0) there is a single critical force /*. 




Figure 5. The temperature dependence of the critical force for horizontal pulling in 
the n oo limit. When the force is less than the critical force the walk is compact 
and when it is greater than the critical force it is expanded. 



At fixed force and below the critical temperature (which depends on the force), we 
have Zc{h, l,uj) = l/u and so 

«:(/,, 0, /?)=/? for P>P\f,,0). (5.3) 

For < P < P^ifx, 0) we know that Zc{h, 1, u) is a strictly decreasing function of P, and 
so K{f^,0,P) is a strictly increasing function of /3. Near the transition we know from 
Owczarek and Prellberg (2007) (see the discussion after equation (3.9)) that 

«:(/., 0, /?) - /5 ~ CiP' - Pf' as /5 ^ ipy (5.4) 

and similarly 

«:(/., 0, /?) - /5 ~ D{f, - fif/' as (/*)+ (5.5) 

for some constants C and D. 



A self-interacting partially directed walk subject to a force 18 

For < (3 < 0), the singularity in z closest to the origin in G{z;h,l,u;) is a 

simple pole. From this we deduce that 

Znifx, 0, (3) = e'^(^-°'^)"+°(^) for < /3 < 0) . (5.6) 

In the semi-continuous model (Owczarek et al. (1993)) the following was derived 
for the scaling of the partition function at low temperatures 

where Ksifx,^, (3) is a negative non-constant analytic function of f3 and fx, and also 
we have that lim^^^t ^^(/a;, 0, /?*) = 0. Both the discrete model (given explicitly by 
Owczarek and Prellberg (2007)) and the semi-continuous model (given by Owczarek 
et al. (1993)) have generating functions that have uniform asymptotics with the same 
algebraic structure, given by a ratio of Airy functions. In particular the results in 
Owczarek et al. (1993) depend on the scaling of the locations of the poles in the 
generating function which has the same scaling in both cases. Without making the 
explicit calculation one can surmise that equation (5.7) also holds for the discrete model. 

For completeness wc note that at /? = (3^{fxj0) there is an algebraic singularity in 
the generating function (sec Owczarek and Prellberg (2007)). Assuming the conditions 
of Darboux's theorem hold gives us 

Zn{f.,0,(3\f,,0)) = e/^'-§Mn)+o(i) _ ^5_g) 

One could calculate the average extension in the generalised ensemble simply by 
differentiating the generating function with respect to the force. In the discrete model 
this gives a complicated expression in terms of g-Bessel functions and its derivatives. 
However as the equivalence of the thermodynamic limits of the canonical ensemble and 
generalised ensemble is limited to < /3 < /3* it is simpler to use the results above to 
deduce the behaviour of the average extension in the canonical ensemble directly. In 
the canonical ensemble the average extension {sx)nifx, fy, P) is defined by 

W.(/.,/.,/3) = i^^^!^^^^. (5.9) 

Assuming the conditions for differentiating the asymptotic expansion term-by-term hold 
the result (5.6) implies that 

{sx)n{fx. 0, (3) = Y^^^Q^'^ n + 0(1) for < /3 < /3* , (5.10) 
so that the thermodynamic limit extension per unit length 

Sxifx, fy: f3) = lim ^'^^r.(fx,fy,P) ^5 -^^^ 

n— >oo n 



IS 



SxifxAP) = ^^^^^ > for < /3 < /3* . (5.12) 
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Moreover the scaling near /3* in equation (5.5) as given by Owczarek and Prellberg 
(2007) implies that 

Sxifx, 0, /3) ^ 0+ as /3 ^ (/3*)- . (5.13) 

One can calculate this quantity from the generalised ensemble in the semi-continuous 
model (see equation (3.46) in Owczarek et al 1993) and analyse it in a more 
straightforward manner than in the discrete case (since it involve Bessel functions rather 
than g-Bessel functions), and it shows the same behaviour. 

On the other hand the result (5.7) (once again assuming differentiability of the 
asymptotic expansion) implies that 

(sMUAP) - l ^^-^^-°'^^ nV^ + 0(l) for/3*</3 (5.14) 

and so that 

5,(/,,0,/3) = for (5.15) 

We have calculated the stress-strain curves by computing the free energy from the 
boundary of convergence (section 4) and then numerically differentiating. Using these 
results and the above arguments we can make the following comments. In all cases the 
function of the average extension per unit length S^ifx, 0, (3) is a continuous function of 
fx for fx > 0. 

At high temperatures where P < /3*(0,0) there is no critical force and the average 
extension per unit length Sx{0,0,P) > at zero force. At such fixed P the function 
Sx{fx,0,P) is an analytic function of fx for fx > and is strictly increasing with 
increasing force. There is a unit horizontal asymptote, approached from below for large 
forces. There is a finite nonzero slope ^'^'^^^j'"-^) at fx — 0. The qualitative behaviour is 
the same as that shown in figure 2 for the the oo = 1 case. 

At the critical temperature for zero force /3 = /9*(0, 0) the average extension per 
unit length Sx{0,0,P) = at zero force. Again the function Sx{fx, 0, P^{0,0)) is an 
analytic function of fx for fx > 0, is strictly increasing with increasing force, and has 
a unit horizontal asymptote, approached from below for large forces. Now however the 
slope ^^^^■f'^'^jf^ (Q''^)) diverges when fx 0+ as fx^^^. This is shown in figure 6. 

At low temperatures where /3 > /9*(0,0) (so there is a critical force /* > 0) there 
are two regimes. For < /j, < /* the average extension per unit length Sx{fx, 0, /3) = 
regardless of the force. However for /* < fx the function Sxifx:^:^) is an analytic 
function of fx, is strictly increasing with increasing force, and has a unit horizontal 
asymptote, approached from below for large forces. Again the slope diverges 
when fx {fi)~^ as {fx — /^)~^^^- This is illustrated in figure 6. 
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Figure 6. The dependence of lim {sx)n/n, the hmiting average horizontal span 

n — >oo 

per unit length, on (3fx for different temperatures: at /3 = /3c (upper curve) and at 
/3 = log 4 > Pc (lower curve). Note that above the critical temperature the behaviour 
is qualitatively the same as that shown in figure 2. 



6. The full model: Pulling in both directions 



6.1. No self-interactions (uj = 1) 

As in subsection 3.1 which describes the solution when there is no vertical pulling (and 
no self-interactions) one can write down a simple functional equation for the generating 
function G(x, ?/_, 1) as 



so that 



G{x,y+,y-, 1) = x(l + 
l + G{x,y+,y-,l) 



y+ 



y- 



1 - 2/+ 1 - y- 
1 



-)(l+G(x,y+,?/_,l)) 



1 — X-, 



i-y+y- 



^ {l-y+)il-y.) 

and setting x = hz, y+ = zv and ?/_ = z/w we get 
l + G{z;h,v,l) 



1 - hz 



l-{v+v-'^)z+z^ 



(6.1) 
(6.2) 

(6.3) 



So the answer is a simple rational function, and the critical z (that is Zc{h, v, 1)) is the 
root of a cubic. For h = v = 1 this simplifies as expected. 
If we set fx = then 

1 + G{z; 1, V, 1) = r \ . T ' (6-4) 



1-z 



1-2^ 



l-{v+v^'^)z+z'-' 
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and from this we can readily calculate the ratio {sy)/{n) and take the thermodynamic 
limit by letting z Zc{l, v, 1). In figure 7 is a plot of lim^^2^(sy)/ (n) against Pfy. 
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Figure 7. A plot of lim {sy) / {n) when ui = 1 against /3/, 



Solution of the full model 

In the appendix of Owczarek et al. (1993) it was shown that if one adds fugacity variables 
I/+ for steps in the positive vertical direction and y- for steps in the negative vertical 
direction then a generalisation of the method for solving the standard problem yields 
the generating function G{x,y^,y_,u) as described in section 2. Using = y^u and 
g_ = y-U it is given by 

1 + iu) = j-^ 1^"!"! TTT , (6.5) 

\21-i{x, y+, y-,uj) - [l + uj + [1 - uj)x)\ 



where 



with 



""^"'^^'^-'"^ ^ iAt + Bt)iAt - Bt) - + B,)iA, - B,) ' (^"^^ 



°° „2»n+l/, , i\2m+l('„ „ Vm{m-|-r)„ r-|-m+l 
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and 

P[A] = (A- l)(A-a;). (6.9) 

Note that q+q- — q^, and given that all the parameters are positive we have q = ^yq+q-. 
This solution is then a clear generalisation of the form for the horizontal pulling model 
where 

1 + G{x,y,y,uj) = J— ^ ^ — — --y. (6.10) 

[2n{y, y(jj, xy{uj - 1)) - (1 + u; + (1 - uj)x)\ 

where TC{y,q,t) — H{y,q,qt)/H{y,q,t) is given in terms of H{y,q,t) defined by 
equation (3.31). 

However the problem of generalising the analysis of Prcllbcrg (1995) is a daunting 
one. On the other hand wc can still make some deductions about the behaviour of the 
generating function. Since the partition function is positive we know that G{z] h, v, ui) 
is a strictly increasing function of z for fixed h,v,u! > 0. 

Let us consider u; > 1. The functions Af and Bf" converge whenever q < 1 and 
moreover there are singularities when \q\ — 1. The solution for the generating function 
G{x,y^,y-,u!) is given in the next section and one finds that G{x,y+,y-,uj) on this 
surface is finite when u > out, where out is given by the solution of equation (6.27). The 
only singularities that can occur in the generating function for g < 1 are poles occurring 
when the denominator of (6.5) is zero. The nature of these poles is unclear but we 
conjecture that they are simple poles and that they only exist when uj < ujf 



6.3. Solving the full model on the special surface z = l/ui 

While expression (6.6) is rather unwieldy, it is possible to find a sub-manifold of 
parameters for which the solution can be made more explicit. If we restrict to 

q+q- — q^ — 1 that is q — cuz — 1 , (6-11) 

then (6.6) has a singular limit. This is analogous to the fact that in the case of (6.10) 
one obtains an algebraic function for g = 1. 

Hence we have q+ = v and g_ = 1/v. Also, we shall retain the variable x for 
convenience during the calculations before substituting x = h/cu as required in the final 
expressions. 

Let g^ = y±~^Gf where Gf is the generating function for walks with r steps in the 
positive and negative directions respectively, and we have in analogy with the derivation 
of the full problem 

g%, + l)(g± + l)g%, + (o;(l + q+^) + {u + lfq±)g%2 (6.12) 

-ujq±{uj + l)(g± + l)g^j^i + uj'^q±^gf = q±x^{uj - lfgf+2 ■ 

The characteristic polynomial is 

P^(A) = (A - 1)(A - uj){\ - g±)(A - q±uj) - q±x\cu - ifX^ . (6.13) 
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Despite being a polynomial of degree four, -P±(A) has sufficient symmetry to allow for 
simple explicit solutions. The key observation is that if A is a root of -P±(A), then so is 
q±uj/X. If we define 

/.± = A+5^, (6.14) 

then ijl± satisfies 

/i| - + l)(g± + l)/i± + {uj + g±)(l + ujq±) - x^{u - l)^ = . (6.15) 

We can therefore obtain all four roots of -P±(A) by solving the quadratic equation (6.15) 
for ij,±, followed by solving the quadratic equation (6.14) for A. 
The general solution for (6.12) is a linear combination 

4 4 

= J2 AiXl,i respectively Gf = J] ■ (6.16) 

i=l i=l 

After some algebra the generating function l+G(l/a;; h, v, cu) — G{h/u!,v/u!, l/vcu, cu) 
can be found to be 

h \ (^1^2^; - l)itiv - l)it2V - l)(a; - l)u;^ 

l + G{l/u;h,v,u) = ^5715 ' ^ 

where 

D{ti, t2) = t-it2{uj - ifh^ + {uj- l){tit2V - l){tit2u;v - l)hw 

- (ti - l)(tiv - l)(t2 - l)(t2V - l)a;3 (6.18) 
with ti and t2 the two distinct roots of L{t) = with 

L{t) = vh'^{(jj - 1)V -{t- l){ujt - l){vt - l){ujvt - \)u? . (6.19) 
We note that L(i) is a quartic polynomial in t. If we let 

s^^^vut (6.20) 
then s satisfies 

- (u; + 1) (^; + l)s + (u; + ^;) (1 + u;^;) = vh^ ^'^~}^ , (6.21) 

and the two solutions of this equation are 

_ (u; + l)(i.- + 1) ± ^{uj \f (^2 _ 2^ + 1 + 4i;/iVcJ^ 

s±— 2 ■ yQ.ll) 

Now only one of the solutions of (6.20) is apphcable, namely 

S — \J — 4:VUJ 

t = ^ 6.23 

2vu; ^ ' 

which gives us that 
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and 

Importantly, on substitution into the generating function the denominator can be found 
to be a polynomial with a factor 

{vuo'^ + 2hu;v + vh"^ - 2vh'^uj + uj^v - 2hu^v + vh'^u'^ - - v'^u^) (6.26) 

and no other relevant factors. Hence the generating function is singular on the curve 
q = 1 when this factor is zero. It can be seen that the algebraic singularity (a square root) 
occurs in the generating function at the same place by considering the discriminant of 
(6.13). It can therefore be deduced that the exponent 7^ = 1/2 describing the singularity 
in the generating function approaching the transition point tangentially (as described 
in Owczarek et al 1993). The exponent is independent of the value of h and v and so of 
whether there is a horizontal and/or vertical pulling force. 
We therefore have the location of the critical point as 

a;^(l + o;^) + 2a;(l - u;^)h + (1 - u;)^h'^ = u;\v + l/v) . (6.27) 

One can rewrite this as 

cosh{pfy/2) = exp{-l3/2) cosh(/3) 

- exp(-/3) sinh(/3/2)(exp(/?/,) - 1) . (6.28) 

When there is no horizontal force, that is /i = 1, we therefore find the critical point to 
be when 

{uj^ + lf = {v + 2+ i)a;3 ^ ^6.29) 

which reduces to the known result of Binder et al. (1990) when v — 1. Prom this 
equation the critical force temperature plot can be found as 

/* = cosh-^ (2e-^ cosh2(/3) - l) . (6.30) 

In figure 8 we plot the critical force against temperature for vertical puUing. Note that 
as the temperature approaches the critical value for no pulling force the slope of this 
curve diverges in contrast to the analogous horizontal pulling curve. 

Note that when v — lin (6.27) we find equation (5.1), that was derived in a previous 
section. 

7. Functional equation method 
7.1. Horizontal pulling 

For the case of horizontal pulling we can construct a recursive functional equation for 
the partial generating function, cf section 2. Let us define 

, . . 2Gr{x,y,y,u) if r > 1 
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1.2n 




Figure 8. The temperature dependence of the critical force for vertical pulling in the 
n — > oo limit. 

SO that 

F{p) = J29rf- (7-2) 

T->0 

The functional equation is constructed by considering what happens when an extra 




Figure 9. Schematic representation of the terms of the functional equation (7.3) 
arising from the five cases described in section 7.1. 

column is added to the left of a walk. Let Vr be the set of all paths with at least one 
horizontal step followed by r > first column vertical steps. Let r' be the number 
of second column vertical steps. We can partition Vr into the following five disjoint 
subsets. 
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Case I: The walk has only one step. This is generated by x. 

Case II: These are all walks with only one horizontal step and at least one first column 
vertical step. These are generated by 2xyp/ (1 — yp), the factor of two giving all upwards 
and all downwards sequences of steps, both generated by yp/ (1 — yp). 

Case III: These are all paths with at least two horizontal steps and the first column 
vertical steps (if any) are in the same direction to the second column vertical steps 
(if any). These are generated by xF{l)/{l — yp). There is no factor of two for the 
following reason. The product xF{l)/{l — yp) corresponds to concatenating a sequence 
of vertical steps, V ^ — yp) to and arbitrary path F(l). The generating variable 
p is set to one in F{p) since the leftmost vertical steps arise from the V sequence. If 
the first vertical sequence of steps of -F(l) are upward, then the V steps are interpreted 
as also being upward (and hence no contacts are created by the concatenation, thus 
no u factor). If the first vertical sequence of steps of F{1) are downward, then the V 
steps are interpreted as also being downward (and again, no contacts are created by the 
concatenation) . 

Case IV: These are all paths with at least two horizontal steps and < r < r' and 
with the further condition that the vertical steps in the first column are in the opposite 
direction to those in the second column. These are generated by xypF {uoyp) / {1 — yp). 
The argument for this form is similar to that of Case III, except now contacts are created 
by the concatenation. The contacts are accounted for as follows. All new contacts occur 
in the overlap between the first and second columns. Since the p in F{p) tracks the 
number of vertical steps in the second column, the new contacts can be accounted for 
by replacing p by ujp. The factor of yp in F{ujyp) is interpreted as giving rise to first 
r' vertical steps in the first column (these can be thought of as 'virtual steps') and the 
remaining r' — r vertical steps are generated by 1/(1 — yp) and are interpreted as being 
concatenated on to the virtual steps, as illustrated by the fourth term in figure 9. 

Case V: These are all paths with at least two horizontal steps and < r' < r and with 
the further condition that vertical steps in the first column are in the opposite direction 
to those in the second column. These are generated by xuiyp [F{1) — F{ujyp)] /(I — ujyp) 
which can be shown using inclusion-exclusion. The term xujyp/ (1 — uyp) counts 
an infinite sequence of vertical bonds and contacts. So the generating function, 
Ti = xu!ypF{l)/{l — cuyp), counts all the required configurations, but also those 
with r' > r. The contribution of these over-counted configurations is then given by 
T2 = xujypF [ujyp) / {1 — ujyp) . Hence the contribution of this case is Ti — T2. See 
figure 10. 
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Figure 10. Schematic representation of contribution of Case V to the functional 
equation. 



Combining all five cases together, as illustrated in figure 9, generates all the walks. 
Thus we get, 

F(p) = x + 2x + X — ^ — F(l) + X F{ujyp) 

1-yp l-yp l-yp 

+ a:-^^^ \f{1) - F{uyp)] . (7.3) 
1 — ujyp L J 

This functional equation can be solved for F{p) by the method of iteration (Bousquet- 
Melou 1996). 

The case uy = 1 is much simpler to solve. Letting uy = 1 gives, 

F{p) = x + 2x + X — ^ — F(l) + X F{p) 

1-yp l-yp 1-yp 

P 

+ Xz 



F(l) - Fip) . (7.4) 
1 — pi J 

Note, to find F{1) from equation (7.4) we cannot put p = 1 because of the denominator 
1 — p. To avoid this problem, first collect coefficients of F{p) to give, 

K{p)F{p) = x{l - p){l + yp) + x{l - yp^)F{l) (7.5) 

where the "kernel", K{p), is given by 

K(p) = yp^ — (\ — X + y + xy) p+1. (7.6) 

If we can put K{p) = in (7.5) then, since the lefthand side is zero, we can solve 
for -F(l). Note, the condition K{p) = implicitly constrains p, which, since K{p) 
is quadratic in p, implies that we get two possible functions p±{x,y). However some 
care must be taken as this way of finding F{p) assumes F{p)K{p) = which is only 
the case if lim^^p^ F{p)K{p) = 0. Thus, we first assume this is the case for one of 
the solutions, say p+ and then verify this assumption once F{p) has been explicitly 
computed. As is readily shown, each of the two assumptions lim^^p^ F{p)K{p) = 
and limp^p_ F{p)K{p) = give rise to different functions, denoted F±{p). In fact, F{p) 
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is an algebraic function and each of the two functions F±(p) are the two branches of 
F[p). The Umit Ump_^p^ F{p)K[p) only vanishes if the correct branch is combined with 
the correct limit (the other branch has a pole at p+ and hence the limit is the non-zero 
residue of F{p) at 

Thus, assuming Mnip^p^ F^{p)K{p) = and taking the limit as p — >• p+ on both 
sides of (7.5) gives F_^_(l). Thus we obtain the branch -F+(l) on the hne ujy — 1 as 

F41) = Jl:iP±}il±P±yl (7.7) 

i-pXv 

or, exphcitly 



l + i^+(l)-W^2(^_l) + 2a;(|/Vl)+y-l- ^^-^^ 

Using this solution, one then verifies the assumption that limp_^p^ F_^[p)K[p) — 0. This 
solution has a square root singularity at 

x\y - 1) + 2x{y + 1) + y - 1 = 0, (7.9) 

which gives the same critical value of x, Xc as given by equation (5.1) using the 
transformations x — h/cu and y = l/cu. Thus we see that, using this method, Xc 
arises via p^ and hence from the kernel. 

The functional equation (7.3) is also closely linked to the recurrence relation (3.16). 
If we make the substitution (2.21) into (7.3) and equate coefficients of p, then, for r > 
we get 

(r-l r-1 \ 

fe>0 fe=0 K>0 k=0 / 

which is readily put in the same form as (3.16). 

7.2. Vertical pulling 
Let us define 

oo 

F+(p) = J2GHx,vy,y/v,u;)f (7.11) 

r=0 

and 

oo 

F-{p) = J^G;{x,y/v,yv,u;)p' . (7.12) 

r=l 

Using similar arguments to the derivation of equation (7.3) we obtain the following 
pair of coupled equations. 

F^ip) = x-\- XV — ^ h X -^^'''(1) + X — ^ F~iujypv) 

1 — ypv i- — yp 1 — ypv 

+ x "^^^ \F-(l)-F-(ujypv)] (7.13) 
1 — ypcuv I J 
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ypv 1 — ypv 
ujypv 



1 — ypv 
F+(l) - F+{ujypv) 



(7.14) 



1 — uypv 

These equations can again be solved by the method of iteration, resulting in complex 
g-series similar to (6.6). 

There is also the simple case occurring when cuy — 1. Using the notation 

L'{a) = , L{a) — aL'{a) and AL{a) — L{ya) — L{a) , 

we can write the two functional equations as 

F+(p) = 1 + L{ypv) + L'{ypv)F^{l) + L{ypv)F~{pv) 

+ L{pv) \f~{1) -F'ipv) 



and 



F (p) = L{ypv) + L'{ypv)F (1) + L{ypv)F^{pv) 
+ L{pv) \f^{1) -F^{pv) 



where 



xF^ 



and 



xF~ 



u>y=l 



ojy=l 



(7.15) 



(7.16) 



Solving the two equations for F (p) and F (p) gives the pair 

K^{p,v)F~^ (p) = 1 + L{ypv) + xL{yp)AL{pv) 

+ [L'{ypv) + xL{p)AL{pv)] F^(l) 

+ [L{pv) + xL'{yp)AL{pv)] F'{1) (7.17) 

and 

K'ip, v)F~{s) = L{ypv) + xL{yp)AL{pv) [1 + L{yp)] 

+ [L'{ypv)+xL{p)AL{pv)\F~{l) 

+ [L{pv) + xL'{yp)AL{pv)] F+(l) , (7.18) 

where the two kernels are given by 

K^{Pi = (1 ~ ~ yp)(l ~ pv){l — ypv) — x'^p^v{y — 1)^ (7-19) 

K-{j),v) = K+(p,v) . (7.20) 
Thus we choose p such that 

K-{p{v),v) ^ K+{p{v),v) (7.21) 

which is the same equation as the characteristic equation (6.13) with p \ and 
y replaced by l/u. Thus we see that the kernel that arises from the functional 
equation approach corresponds to the characteristic equation (6.13) required to solve 
the Temperley recurrence relations (6.12). 
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8. Three-dimensional model 

The model can be generalized to a three-dimensional partially directed walk model in 
which the walk is self-avoiding, can take steps in the +x or +y directions or in the ±z 
directions. Again it will be convenient to require that the first step is in the -\-x or -\-y 
direction. We shall only consider a force applied in the {x, y) -plane and as such we shall 
keep track of the span in the x and y directions as well as the number of contacts. 
If we write c„ for the number of these walks with n steps it is easy to see that 

lim n-^ log On = log[(3 + VT?) /2] . (8.1) 

n— >oo 

Suppose that Cn{sx, Sy,m) is the number of walks with n steps, m contacts, and with 
spans Sx and Sy in the x and y directions. The corresponding canonical partition function 
is 

Z,,{K,hy,uj) = J2 Cnis^,Sy,m)hl^hlyu"' (8.2) 

and we define the generating function 

G[z\ hx, hy, (^) = ^ Zn{hx, hy, uj)z^. (8.3) 

n 

Concatenation arguments can be used, as in section 3, to establish the existence of 
the limiting free energy 

K{h^,hy,uj)= lim log Zn{h^,hy,uj). (8.4) 

n— >oo 

Methods exactly analogous to those in section 3 can be used to show that ^{hj;, hy,uj) 
is convex and continuous, and differentiable almost everywhere. At fixed uj, and hy 
G converges if 2; < Zc{hx, hy,Lj) = exp[—K,{hx, hy,Lj)], and the phase boundary z = Zc is 
continuous. 

If we turn off the interactions by setting cu — 1 and write Gq = G{z;hx,hy,l) then 
Go satisfies the equation 

Go = {hx + hy){Go + l)z + ^"''^/'-^^''^^ Go + 1) (8.5) 

i. — z 

so that 

I — z 

^"'^'^ i-z-z{z + i){K + h,y ^^-^^ 

The force-extension curve can be calculated as in section 3 and this has the same general 
form as for the two-dimensional model, and as found experimentally for good solvent 
conditions (Gunari et al. 2007). 

For attractive interactions (a; > 1) the Temper ley approach described in section 3 
for the two-dimensional case immediately generalises to this model. Defining partial 
generating functions as in section 2, the partial generating functions obey the relations 

go^{K + hy)z{l + G) (8.7) 
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and 

gr = {h., + hy)z''+' (2 + J2i^ + ^')9k + {l+u;nJ29k] ^ (8-8) 

\ k=0 k>r / 

SO that Qr satisfies the relation 

gr+i - + q)gr + q''iK + hy)^iq - ^)9r + qzgr-i = o (8.9) 

with q = ujz. 

The analysis goes through as in section 3 and one can find the temperature 
dependence of the critical force. By choosing the relative magnitudes of and hy one 
can change the direction in the (a;, ?/)-plane in which the force is applied. In all cases the 
force-temperature curve is monotonic and does not show reentrance. Presumably this 
is because the model does not have extensive ground state entropy. The phase diagram 
can also be calculated, using the same methods as in section 3, with qualitatively similar 
results. 

1.2 n 



1.0- 



0.8- 



0.4- 



0.2- 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 

Figure 11. The temperature dependence of the critical force for the three-dimensional 
model in the n — > oo limit when the force is applied in the a:-direction. 

If we set /ly = we reproduce the results of section 3. We consider two other 
cases. First we set = h and hy = 1. This corresponds to applying a force in the 
x-direction. We show the critical force-temperature curve in figure 11 and the boundary 
of convergence in figure 12. 

These are qualitatively similar to the corresponding figures for horizontal pulling 
in two dimensions, though with quantitative differences, of course. In particular the 
critical force-temperature curve has zero slope in the T — > limit and is monotone 
decreasing. We also set = hy = h which corresponds to pulling in the xy-plsme but 
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-I • 1 • 1 ■ 1 • 1 • 1 • 1 • 1 

1 2 3 4 5 6 7 

UJ 

Figure 12. The boundary of convergence of the generating function as a function of 
w for the three-dimensional model when the force is apphed in the x-direction. The 
rectangular hyperbola is independent of the value of h^. The three curves with points 
marked correspond to hx = I (top curve), = 1-2 and hx — 1-5. 
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0.8 



0.4 
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Figure 13. The temperature dependence of the critical force for the three-dimensional 
model in the n oo limit when the force is applied in the xy-plane at 45° to the x-axis. 

at 45° to the x-axis. We show the critical force-temperature curve in figure 13 and the 
boundary of convergence in figure 14. 

The critical force now has negative limiting slope in the T limit but remains 
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Figure 14. The boundary of convergence of the generating function as a function of 
Lu for the three-dimensional model when the force is applied in the xy-plane at 45° to 
the a;-axis. The rectangular hyperbola is independent of the value of hx = hy. The 
three curves with points marked correspond to hx = hy = 1 (top curve), h^ = hy = 1.2 
and hx — hy = 1.5. 



monotone decreasing. This is because of the curious feature of this model that the 
ground state has no (extensive) entropy in the compact state but acquires entropy 
under the influence of a force. The low temperature behaviour can be understood 
by the following crude low temperature argument. Think of an n-edge walk at low 
temperature T under a tensile force / at 45° to the x-axis. li n — m edges of the walk 
are in a compact state and the remaining m edges are extended the (extensive) free 
energy can be written as 



where e < is the vertex-vertex attractive energy. Differentiating with respect to m 
and setting the derivative equal to zero gives 



Setting e = — 1 gives a critical force of 1 at T = and the limiting slope d//dT = 
-log2 < 0. 

9. Discussion 

We have analysed the polymer model of partially directed walks with self-interaction, 
so as to induce a collapse transition, under the influence of tensile forces on the ends of 
the polymer. The problem of forces only in the preferred direction of the walk can be 





/ 



e -Tlog2. 
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elucidated completely. The phase transition which is a second-order transition without 
any force is unchanged by the presence of such a force. The force extension curves at 
high temperatures look qualitatively similar to those of AFM experiments (Gunari et 
al. (2007)). 

The solution of the full model is more problematic. While the exact solution of 
the generating function of partition functions can be written down in terms of g-series, 
these functions are even more complicated than those encountered in the standard model 
that has only forces in the preferred direction (those are gr-Bessel functions). We have 
been able to solve the model on the important special curve in parameter space which 
should contain the transition point. This seems to indicate that, again, the transition is 
unaffected by the force. This is more difficult to understand physically as the force must 
change the high temperature state of the polymer, unlike a force applied only in the 
preferred direction. It would therefore be interesting to analyse this full model further. 
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